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We use light-front dynamics to calculate the electromagnetic form-factor for the Hulthen model of the deuteron. 
For small momentum transfer Q 2 < 5 GeV 2 the relativistic effects are quite small. For Q 2 ~ 11 GeV 2 there is 
~ 13% discrepancy between the relativistic and non-relativistic approaches. For asymptotically large momentum 
transfer, however, the light-front form factor, ~ log Q 2 /Q 4 , markedly differs from the non-relativistic version, 
~ 1/Q 4 . This behavior is also present for any wave function, such as those obtained from realistic potential 
^5 ' models, which can be represented as a sum of Yukawa functions. Furthermore, the asymptotic behavior is 

in disagreement with the Drell-Yan-West relation. We investigate precisely how to determine the asymptotic 
behavior and confront the problem underlying troublesome form factors on the light front. 

o 

1 Introduction 
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The light front approach to quantum dynamics was introduced by Dirac jl]] a half century ago. Since then, light front 
dynamics has developed into an active area of research for a variety of reasons, e.g. its minimal set of dynamical 



operators, the simplicity of the light-front vacuum, and the close connection to experimental observables. Light front 
techniques have long been used in analyzing high energy experiments with nuclear and nucleon targets ||, [|, ^, ||] . 
Indeed light front dynamics is relevant to describe such reactions, since, for example: in the parton model, the ratio 
k + / p + (where k + = fc° + k 3 is the plus momentum of the struck quark and p + that of the target) is nothing more 
than the Bjorken x- variable. 

Some recent efforts have been made || to render the theory more understandable by using models reminiscent 
of basic quantum mechanics rather than by invoking quantum field theory. These models find particular reality 
in nuclear physics 0, where some nucleon interactions may be described by a mean field potential. Nonetheless, 
the similarity of the light-front bound state equation to the Schrodinger equation is grounds enough to put familiar 
quantum mechanical problems on the light front. Below, we do precisely this for the Hulthen model of the deuteron 
and its electromagnetic form factor. Of particular concern here is the asymptotic behavior of the form factor, which 
differs from the non-relativistic version. This may be of interest to experimentalists seeking to probe asymptopia. 
Recent measurements of deuteron form factors at the Jefferson National Laboratory || reach momentum transfers 
of Q 2 = 6 GeV 2 and future projects hope to reach upwards of Q 2 — 11 GeV 2 . In this range of momentum transfer, 
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there is ~ 13% discrepancy between the relativistic and non-relativistic form factors calculated in this paper (as we 
will illustrate in figure |^) . 

This paper's organization is similar to that of a detective story. First we recall a minimal amount of light 
front dynamics in section and explain how we put the non-relativistic Hulthen potential on the light front. Next 
in section we calculate the electromagnetic form factor using light front dynamics and compare with the non- 
relativistic version calculated in section ^. The low momentum behavior of the these form factors shows only minimal 
differences, while the high momentum behavior leads to surprising trouble in asymptopia (section |5|). We could solve 
the mystery at this point by deriving the asymptotic behavior of the form factor. Instead, we proceed by assuming 
factorization holds in the asymptotic limit. This leads us to consider various previous attempts at dealing with the 
end-point region and to dispel any lingering misconceptions. In section |^, we discover that troublesome asymptotic 
behavior also lurks in other models on the light front. With enough clues at hand, we are able to pin-point the cause. 
The asymptotic behavior is then deduced in section 0, and is similar to that obtained for the Wick-Cutkosky model 
in Ref . || . Finally, we summarize our findings in a brief concluding section. 



2 Hulthen model on the light front 

In light front dynamics, one quantizes the fields at equal light-front time specified by x + = x° + x 3 = t + z. This 
redefinition of the time variable leaves us with a new spatial variable x~ — x° — x 3 — t — z. The remaining spatial 
variables are left unchanged by this transformation: = (a; 1 ,^ 2 ). 

If one uses x~ as a spatial variable, then its momentum conjugate is p + = p a + p 3 . This leaves p" = p° — p 3 as 
the energy, or the x + -development operator. The details of this formalism do not concern us here — the interested 
reader should consult [jl0| for a good overview. What is important to note, however, is that the relativistic dispersion 
relation p^p^ — m 2 takes the form 

2 o 

p = — ^ w 

and thus the expression for the kinetic energy avoids the historically problematic square root. 

For a bound state of two particles interacting via a potential V, the light-front wave function is determined by 
solving the equation 

^ = : ~ 1 fcXCT ^' (2) 

where M is the invariant mass of the system, m, the particle mass, and xi is the plus momentum fraction carried by 
the i th particle, namely Xi = kf / P + , with P + as the total plus momentum, k^ +kj. Let us take the nucleons to be 
of equal mass, and use m — T "p+ m " as the nucleon mass. Furthermore, since we have only two particles, the sum of 
x\ and X2 is one. So we choose x\ = x and consequently, x-i = 1 — x. 

In order to simplify Eq. ^|, it is customary to define the relative light-front variables ]T^ | 

P ± =k^+k 2 L (3) 
p 1 - = -xkj L +(l-x)k^-. 



Straightforward algebra transforms Eq. ^| into 

2 ,A = ( 

x{\ 



M 2 ^=( P± n +? " 2 +V)r. (4) 
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which is the coordinate representation of the Weinberg equation jl3|]. Eq. |] is still quite complicated to solve, so we 
define an auxiliary operator 



p 3 



to cast the equation into a familiar form. Defining M = 2m — e (where e is the binding energy) and using the above 
definition, wc find 

(j - em)v = (p^ 2 +(p 3 ) 2 + j)v> 

^( P 2 + v H )^ (6) 

where we have efficaciously chosen V H = V/4 to be the Hulthen potential [ O . 

The above equation is the coordinate-space version considered by others, see e.g. jn|. Taking p conjugate to r, 
we have 



and the well known ground state solution 

ip(r) cx (8) 



e -ar _ £ -br 



with a = \J em — e 2 /4 as dictated by Eq. |^. The experimentally determined values of the model parameters are [ [T(i[ : 
a = 0.23161 fm" 1 and b = 1.3802 fm" 1 . 



3 Electromagnetic Form Factor 

The electromagnetic form factor on the light front has the form [[l7) 

2 r dxdp^ , ± x pX 

J x(l-x) 

where the momentum transfer Q 2 = q^ 2 . The momentum-space Hulthen wave function is the Fourier transform of 
our solution Eq. ||, namely 



mVN / 1 1 



4 V a 2 + p- 1 - +{p 3 ) 2 b 2 + p J - +(p 



■ i \2 



mV Nx(l — x) 



4x(l - x)a 2 + (2x - l) 2 m 2 + p ±2 



(io) 



To calculate the form factor, we must perform three integrals. Writing dp 1 = p 1 dp 1 d<fi with <fi as the angle 
between and q 1 , we see the <j> integral and subsequently the p 1 integral can be computed analytically. Performing 
these integrals leaves us with 

F{Q 2 ) = J f(x,(l-x) 2 Q 2 )dx, (11) 
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Figure 1: Plot of the light front Hulthcn form factor Eq. |TT] as a function of Q in GeV. 



where 



with 



g a /3(x,k 2 ) 



f(x, fc 2 ) = m 2 Nx(l - x) g a0 (x, fc 2 )(« - 6*8? - + S?5g) (12) 
fc 4 + fc 2 (2 j a (x) +Jp(x) + p a p{x, fc 2 )) +j p (x)(p al3 (x,k 2 ) - jp(x) - -y a {x)) 



p a p{x, k 2 ) 
where we have defined 



log 



7 Q (x)( - k 2 + p a p(x, fc 2 ) +7^j(x) - 7q(x)) 

7^(x) = Ax(l - x)p 2 + (2x - l) 2 m 2 
p al 3(x,k 2 ) = \Jk A + 2k 2 jp(x) + 7a(x) 2 . 

We choose N so that the form factor is normalized, F(0) = 1. The constant N is determined by setting Q = in 
Eq. Bl which yields N — 14.931 -1 . Figure [l] shows the form factor as a function of Q. We have also calculated the 
derivative of F(Q 2 ) in the limit Q — > in order to find the root-mean-square deuteron radius 



/ I 



/ dF 
lim = 1.9467 fm. 

Q*-o V dQ 2 



(13) 



4 The Non-Relativistic Limit and Comparison 

Our solution to the Hulthen model on the light front closely resembles the non-relativistic treatment. In fact, we 
have used the non-relativistic solution as a guide in constructing the relativistic wave function. Clearly relativistic 
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effects are contained in the light-front variable x. Quite simply then, the non-relativistic limit of Eq. || is found in 
the limit m — > oo by retaining terms to C[^]- Inverting Eq. |] yields 



-(p 3 ) 2 + m 2 



2m' 



(14) 



Since the measure dx dp 1 ^ — > dp /2m is already first order, we need only keep leading order terms in the wave 
functions to find the non-relativistic form factor. It is clear that to leading order: tp(x, p 1 -) — > ^(p), where the latter 
is the non-relativistic wave function. Quite similarly, we see 



*P(x,p ± +(l-x) q^) 



1 



a 2 + p 2 + p^ • q 1 - + q- 1 - /4 
Returning to the expression for the form factor, the non-relativistic limit is then 



F(Q 2 ) -> / dp^ip + q 1 - /2)^(p). 



(15) 



The form factor above depends neither on the orientation of q^, nor of p. Let us then rotate our coordinate system 
so that q^ is no longer completely transverse. This three-dimensional rotation is only possible now because we are 
integrating dp, which can not be done in the light-front version. Thus we have sent q^ — * q while maintaining the 
length, q 2 = Q 2 . After this rotation, the form factor is strikingly non-relativistic and can be computed analytically 



using Eq. || 



F(Q 2 ) NR = J dr |V(r)|V lq - r / 2 



mN' 



tan 



-i (Q 



(-) 

\4aJ 



2 tan -1 



2(a + b) 



tan 



-x(Q 



(-) 



(16) 



with N' chosen to make F(0) NR = 1. 

From this analytical result, the rms. radius can be easily calculated (see Eq. 



Urn 



Va 4 + 5a 3 b + \2a 2 b 2 + 5ab 3 + 6 4 
2 3 / 2 a6(a + b) 



1.9395 fm. 



(17) 



Comparing with our previous result, the relativistic system is larger by only 0.37%. This confirms our suspicion that 
relativistic effects in this deuteron model are small. We can further confirm this by looking at the difference of the 
relativistic form factor Eq. |ll| and the non-relativistic version Eq. [l6| The percent difference is plotted for low Q in 
Figure || illustrating a difference only ~ 1% in this momentum regime. The small nature of relativistic effects was 
noted early on pS| . 



5 Trouble in Asymptopia 

The above plot shows the absolute percent difference continually increasing as Q increases. In this section, we 
investigate how the light-front form factor compares with the non-relativistic version for large Q. 
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Figure 2: Relativistic and non-relativistic form factor comparison: percent difference: 100 x (F NR • — F)/F NR plotted 
as a function of Q in GeV. 



5.1 Exploring Asymptopia 

Given our analytic expression for the non-relativistic form factor Eq. [H| it is simple to Taylor expand about Q = oo 
to find its asymptotic behavior. To leading order 

lim F(Q 2 ) ~ Mab ^ = 0-080585 (GeV)^ 

The asymptotic behavior of the relativistic form factor is found with the aid of the Drell-Yan-West relation 
(under the assumption that the end-point region dominates the form factor for large Q). This relation takes the form 

lim x f(x, 0) ~ (x - I) 26 ' 1 <^> lim F(Q 2 ) ~ (Q 2 )" 5 . (19) 

The x-distribution function f(x,0) can be calculated analytically using Eq. ^ with = and subsequently 
expanded about x = 1. The leading-order term in the expansion is 0[(x — 1) ], from which we deduce X/Q behavior 
in asymptopia. Given that there were only small differences between the relativistic and non-relativistic form factors 
for low Q, we might expect agreement in the asymptotic region. Moreover, both form factors go as l/Q 4 for large 
Q. So when scaled by Q 4 , at worst the form factors will tend to some common difference as Q — > oo. 

To compare the asymptotic behavior, we have plotted the relativistic and non-relativistic form factors (scaled 
by Q 4 ) for large Q (in figure ^ we plot for experimentally relevant Q 2 , whereas in figure |] we are mathematically 
contrasting the asymptotics). The non-relativistic form factor lines up well with the asymptotic limit predicted by 
Eq. [l^. The relativistic form factor, however, markedly differs from its non-relativistic counterpart, in disagreement 
with the Drell-Yan-West relation Eq. Ft9l. We remind the reader that the relativistic form factor is computed exactly 
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Figure 3: Relativistic and non-relativistic comparison for experimentally relevant Q 2 : the form factors are scaled by 
Q 4 and plotted as a function of Q 2 in GeV 2 . 
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Figure 4: Relativistic and non-relativistic comparison in the asymptotic limit: the form factors are scaled by Q 4 and 
plotted as a function of Q in GeV. 
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for our model [ pp[ . The huge disparity between the non-relativistic and relativistic results, shown in figure |] warrants 
a complete journey through asymptopia. 

5.2 Sharp Peaks at x = 1 

Before proceeding, we note that our light front wave function Eq. [l^ is properly behaved: 

p ± "^(x,p ± ) — > as p ±2 — » oo and 

ip{x,p x ) ->0 as x -> 0,1 (20) 

These conditions stem from non-relativistic versions, and will be trivially satisfied for light-front wave functions 
created using Eq. | and §. Thus knowing the non-relativistic wave functions are peaked for small momenta, our 
light- front Hulthcn wave function must be peaked for small transverse momenta. For large momentum transfer Q, 
Eq. |^ shows that large momentum flows through either ip or ip*. Following Brodsky and Lepage [pT[ , the dominant 
contributions to the form factor in the asymptotic limit come from the two regions which minimize wave function 
suppression 

i) | p ± | <C |(1 — x) q x |, where ^/'*(p" L +(1 — x) (\ x ) is small and V^P" 1 ) is large 

ii) I P 1 " +(1 - x ) q X I < 1(1 - x) q x |, where ^(p^) is small and f/'*(p ± +(1 - ^) q^) is large (21) 

Working first in region i, we can neglect p 1 - relative to q 1 - in ip* since the light-front wave functions are peaked fol- 
low transverse momenta. The contribution from region i is exactly the same as from ii which is made obvious by 
shifting p x . Thus dominant contributions to the form factor in the asymptotic regime appear as 

F(Q 2 )^2 J -^—^*( x ,(l~x)^) J dp 1 - ^(x,p x ) (22) 

4x(l - x)b 2 + (2x - l) 2 m 



8nNm 2 (b 2 -a 2 ) f 1 dx x 2 , 
(1 — x) z 



4x(l - x)a? + (2x - l) 2 m 2 



1 dx g(x) . . 

-ir> (23) 



where we have retained the same normalization constant that appears in Eq. [l2|. [] But to determine the asymptotic 
behavior, we must perform the integral over x which diverges! The end-point region is too peaked for the integral to 
converge, as illustrated by Figure |[ 

The end-point region appears to dominate the form factor in asymptopia (as suggested in [pl[ ). Perhaps the 
actual asymptotic behavior can be extracted by putting the end-point region under scrutiny. Above, we merely 
assumed the validity of the Drell-Yan-West relation; now, we will rigorously investigate it for our model. 

In ascertaining the dominant contributions to the form factor in asymptopia, we have neglected the case x = 1. 
The form factor includes contributions from the end point, but this is where the scheme set up in Eq. ^l] breaks 
down. So the above approximation Eq. |2^ is really only valid for x < 1 — A^, where A is some dimensionless cutoff 
less then one. For 1 — A^ < x < 1, we must return to the full expression for the form factor to get the end-point 

contribution. To leading order, however, (1 — x) q~*~ ~ Amq 1 /Q in the end-point region and the contribution to the 
form factor reads 

F(Q 2 ) EP « / — — f dp x ftrp 1 ^ 1 )^^ 1 ), (24) 
Ji-A-g. x(L - x) j 

1 We obtained the same result by brute force Taylor expansion of Eq. [l^ about Q = oo. Since the integral over x diverges, the series 
expansion of f(x,Q 2 ) lacks uniform convergence. 
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Figure 5: Plot of g(x) appearing in Eq. |2^. The singularity at x = 1 is too severe to bound a finite area under the 
curve. 



where q = /Q is the direction of q^. Since this contribution to the form factor depends only on Q 2 , we can 
rotate our coordinates about the z-axis to make q 1 - parallel to x. The resulting functional form is entirely similar to 
the full form factor, enabling swift evaluation: 



2\EP 



1 



f(x, A 2 m 2 ) dx 



(25) 



with f(x,k 2 ) given by Eq. 12. The leading-order contribution to the above integral in asymptopia is found by 
expanding the integrand about x = 1 and integrating. The result is 



F(Q ) = 8tt7V(6 2 - a 2 )' 



A 6 + 2A 4 - 8A 2 - 2A(A 2 + 1)VA 2 + 4 



log 



Vx 2 



4- A 



A 3 + 3A+(A 2 + l)VA2+4 



(26) 



Q 4 (A 2 + 1)3 

which agrees with the Drcll-Yan-West relation. To quantitatively consider the contribution from the end-point region, 
we have plotted the percent contribution to the form factor from < a; < 1 — A^. We have chosen the value of A to 
be smaller than one. Figure ^| shows that the bulk of the form factor does not come from the end-point region. 



5.3 Divide and Conquer 

In the process of trying to deduce the asymptotic behavior, our model has become infinitely sensitive to the end-point 
region. From the exact Hulthcn form factor Eq. O, however, we know the end-point region does not overwhelm- 
ingly contribute (as Figure || confirms). Our dilemma sounds familiar [p2[ and our approach, not surprisingly, is a 
regularization cutoff. 

To start, let us just toss away the troublesome divergent part of Eq. |2^ by introducing the cutoff A into the x 
integral 



dx 
x(l — x 



-r(x,(l-x)q ± ) / dp X ^fop- 1 ) 



(27) 
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Figure 6: Percent contribution to the form factor (Eq. |ll]) from the region < x < 1 — A^J as a function of Q in 
GeV for different A. As indicated, the end-point region is not dominant for asymptotic Q. 



The cutoff integral above can be computed analytically. Using the integrand of Eq. E3L the result reads 



lim F(Q 2 ) 



32nN(b 2 - a 2 ) 2 

Q 3 



X(a,b,m) 



log A + log — 

m 



(28) 



where x( a ib, m) is a rather complicated, page-long function independent of A. For our parameters a, b and m, we 
have x = 1-3635. Nonetheless, we have discovered behavior via regularization which differs from 1/Q 4 in a model 
which knows nothing about ultraviolet divergences, renormalization, etc. 

One would think that with Eq. |2^, we have determined the asymptotic behavior of the form factor. Although 
we threw away the end-point region to arrive at the above expression, we know precisely its contribution for a given 
A, cf Eq. |2^. The question remains: have we found all contributions to 0[1/Q 4 ]? 

There are order 1/Q 6 corrections to the integrand of equation 23. Putting these in give a correction term: 



-167rm 2 AT(& 2 - a 2 ) 

Q 6 



1-M 



dx x 



(l-xy 



2x(l - x)(a 2 + b 2 ) + (2x - l) 2 m 2 



log 



ix(l - x)b 2 + (2x - l) 2 m 2 
4x(l - x)a 2 + (2x - l) 2 m 2 



Evaluating this correction term to leading order gives: OZ7rJ Q^ 2 - 1 . Thus terms in the integrand of order 1/Q 6 
give a contribution of order 1/Q 4 to the asymptotic form factor. We haven't exhausted all of the 1/Q 6 corrections, 
however — we originally took only the first term in the Taylor expansion of Eq. ^ about p- 1 = (1 — x) q -1 and the 
next non- vanishing term gives contributions of order 1/Q 6 . Even if we were to collect all the C[l/Q 6 ] corrections 
to the a;- integrand, we would have only just begun. One can easily find terms in the integrand of order 1/Q 8 which 
emerge from the regularized x- integral 1/Q 4 . In fact, the integrand's correction terms of any order contribute to 
leading order in asymptopia! 
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Certainly we cannot hope to evaluate infinitely many leading-order terms. At least we have stumbled on to a 
prediction for the functional form in the asymptotic limit. Namely, we have seen 



lim z 4 F(z 2 ) = a + (3 log z, (29) 

Z — >OG 

with z = Q/m. We can test this prediction against the actual form factor's asymptotic limit calculated from equation 
O. In figure |?], we test this hypothesis for empirically determined coefficients a — 0.039472 and j3 = 0.044930 
(calculated for z around 400). As the figure shows, this is indeed the form factor's behavior in asymptopia. It is 
quite curious to note: using Eq. |2^, we would predict f3 = 0.045020 a difference of only 0.20% when compared with 
the empirical value. We believe this discrepancy results from approximating asymptotic z as around 400, not from 
ignoring infinitely many leading-order corrections. Indeed, we never found corrections of order log z/z 4 above, only 
a myriad of l/z 4 terms. It is our belief that the coefficient /3 can be ascertained from the regularization integral 

51 



Eq. E8L The leading correctionn to f3 determined graphically is C[l/z 5 ] — which gives a relative correction of 0.25% 



'0 

for z ~ 400. Taking z larger in order to reduce this term only results in appreciable error in numerical integration. 
To verify our conjecture, we have attempted to find the coefficient f3 by varying z . Figure |8| shows a plot of the 
graphically found value of /3 as a function of z (the midpoint of our interval). Specifically we use a simple linear fit 
in the plot 

(3{z) = \{z + 50) 4 F(z + 50) -(z- 50) 4 F(z - 50) ) log . (30) 

V / z + 50 

The plot shows our cited value (3 = 0.044930 at z = 400. The trend is clear, (3 is increasing to some limiting value 
as z increases. The numerical integration, however, becomes unreliable to ~ 1% past 600. Nonetheless, it appears 
we can determine (i from the regularization integral ^8|. We are at a loss, however, to predict a: there are simply an 
infinite number of correction terms to 0[l/<3 4 ] to evaluate. 

6 Similar Problems 

The problems encountered above are not unique to the Huthen model. In this section, we begin by exploring another 
model with similar behavior in asymptopia. 

6.1 Coulomb Potential 

Let us suppose our particles interact via a Coulomb potential V{r) ~ 1/r, for which we take ip{v) = e~ Kr . Then our 
momentum-space solution to Eq. O is given by 



± m 3 VN 
ip{x,p ) 



16(k 2 + P 2 ) 2 

m 3 VNx 2 (l - x) 2 
(4x(l - x)k 2 + p 4 - 2 +{2x ~ l) 2 ™ 2 ) 2 ' 



2 Expanding the integrand of the form-factor in powers of Q, there are only even terms. Once we exit the cutoff integral, however, we 
can now have any power of Q in the expansion. 
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Figure 7: Asymptotic Behavior of F(Q 2 ) given by Equation 0: W = z jgfejr 9 is plotted as a function of z = Q for 
the parameters a = 0.039472 and (3 — 0.044930 which were graphically determined for z around 400. 
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Figure 8: Simple numerical determination of (via Equation |3Tj) plotted as a function z = Q. The numerical 
integration clearly becomes imprecise past z — 600. 
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where we have used Eq. |5| to re-express p 3 in the light-front center of mass. To find the asymptotic behavior of this 
model's form factor, we use Eq. E3. Performing the integration, leaves us with a logarithmically divergent x-integral: 



lim F(Q 2 ) 



m 6 Nn 



dx x 



Q 4 J {l-x)[4x(l-x)n 2 + (2x-l) 2 m 2 }' 



(32) 



Restricting x to the range < x < 1 — A|j, and performing the integral yields 



lim F(Q 2 

Q— >oo 



m 4 NTr 



x(k, to) — log A + log — 

TO 



(33) 



6.2 Why Regularization? 

As we have seen above, determining the exact asymptotic behavior of light-front form factors is no trivial task. 
The Drell-Yan-West relation is still apt at describing the contribution from the end-point region, however, this 
region does not dominate our form factors for asymptotic Q. Furthermore, techniques to determine the asymptotic 
behavior (brute force Taylor expansion, finding contributions from regions of minimal wave function suppression) led 
to logarithmically divergent ^-integrals, suggesting the non-commutativity of the limits: Q — > oo and x — ► 1. Here 
we show how the potentials we use cause the peculiarities for x — ► 1. 

Starting with the expression for the light-front form factor Eq. ^, we were led to the dominant contribution in 
asymptopia via isolating the regions of minimal wave function suppression, namely 

lim F{Q 2 )^2 f , dX ^(x,(l-x)q ± ) [dp x ^(x,p ± ). 

Q^oo J X{l-X) J 

Now let us utilize the Weinberg equation |l3| , which is the momentum-space version of our Eq. ^: 

^p X ) = Ar m2/ X m > J ^^vMm^yM), (34) 

x{l— x) 

with V as the Fourier transform of the potential. We can use this information in the asymptotic limit of the form 
factor, namely for 

m 2 x f dy dk 



i/>(x, (1 - x) q^) « —2 / -f -My, ^)V(x, (1 - x) q^; y, k^) (35) 

27r d (l - x)Q z J ?/(l - y) 

Of course, in asymptopia the integral containing V can be simplified. The potentials considered above are of the 
form V = V( |p — k| ), where 



2 



,2 



+ m l 



(p-k)-^-k-j +^-V2)y^ y -(y-i/2 ) y^ r j (36) 

which makes explicit use of having two equally massive particles in the center of mass frame ( cf Eq. |5|). As a result of 
the above equation, V(x, (1 — x) q^; y, k ± ) = V(x, (1 — x) q x ; 1/2, 0) to leading order. Henceforth we shall abbreviate 
V{x,{l -x)q ± ; 1/2,0) = V( |q(x)| ), where 



Q 



\q(x)\ = ^{l-x) 2 Q 2 + {q 3 ) 2 = I A /-— (37) 
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Revising the expression for the asymptotic form factor, we find (neglecting overall constants) 



<& nQ^-jj^yvU^), (38) 

where <f)[x) = J dp- 1 ip(x, p- 1 ) and the y-dependence integrated itself away. Eq. ^ contains the answer to our 
troubled journey through asymptopia. At first glance, the integrand appears singular at the end point: containing 
one factor of 1 — x from the measure, and another from the (1 — x) contained in the form factor. These factors 
are quite general and contain nothing specific about the interaction. While Eq. ^ spells out the criteria for good 
wave functions, it is necessary to be further restrictive by requiring ip(x, p^) ~ (1 — x) 2 as x — + 1 if we wish to cancel 
the potentially singular denominator in Eq. pS. This isn't much of an imposition; both the Hulthen wave function 
Eq. and the Coulomb wave function Eq. |3l[go like (1 — x) 2 as x — > l! 

Only one x-dependent piece of Eq. [58] remains to be considered — the potential. It is now immediately obvious 
that the Coulombic form factor should suffer logarithmic divergences in asymptopia: 

1 At 

The potential brings along the anticipated factor of Q 2 , but on the light front, an unwanted (1 — x) tags along. Given 
the behavior of the light front wave function, this extra factor is just enough to make the integral in Eq. |38| diverge. 
The same is true for the Hulthen potential Eq. ^ since we have 

.. t -h/, , m> 12n{a + b) 48n(a + b)x . An . 

Jim yH ( Hx)l ) = - - = - {1 _ x)Q2 ■ (40) 

As illustrated above (Figures 1-4), the form factor itself isn't singular at asymptotic Q, just our means of obtaining 
it. This is made obvious by commuting the limits. Above we looked at Q — > oo first and found problematic behavior 
for x — > 1 stemming from the potential. On the other hand, consider taking x near one first. We already did this 
for the Hulthen form factor, see Eq. Now we are in a position to generalize this result. Since we know our wave 
functions tp(x,p ± ) ~ (1 — x) 2 as x — > 1, the contribution from the end point becomes: 

F(Q 2 ) EP - f {l-xfdx^^-. 

This is just the Drell-Yan-West relation, which, as we have seen, does not account for the majority of the form 
factor in asymptopia. Clearly different behavior is seen when looking at asymptotic Q near the end point, versus the 
end-point region for asymptotic Q\ 

Logarithmically divergent form factors in asymptopia need not plague us any longer. The culprit has been 
unmasked: potentials in the auxiliary coordinate-space, such as the 1/r of the Coulomb model or Eq. ^ for the 
Hulthen, which are ~ 1/r for small r will lead to logarithmic divergences in the expression for the asymptotic form 
factor. Of course, the asymptotic form factor itself is not singular. The logarithmic divergence of our asymptotic 
expression is a thorn-like warning: the series expansion in 1/Q does not converge uniformly in x. 

6.3 Realistic Models on the Light Front 

Above we have seen rather simplistic models lead to electromagnetic form factors with non-standard asymptotic 
behavior. It is not likely, however, that this behavior is physical — though certainly it is the true asymptotic behavior 
for the models considered. 
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Early work on factorization in quantum chromodynamics showed similar logarithms appearing in asymptopia |23] 
for the nucleon electric form factor, and thus a failure of renormalization group techniques. Quite soon there after, 
it was realized pij that these logarithms were just manifestations of neglected higher-order corrections. Apparent 
end-point singularities are removed when evolution of the longitudinal momentum amplitude (f> is properly included, 
and consequently factorization is saved p5| . 

Obviously our models do not have such higher-order corrections, and thus the logarithm remains. We must then 
wonder whether factorization breaks down for more realistic models on the light front. Let us then consider a more 
sophisticated model of the deuteron arising from meson- theoretic potentials |p6| . The general parameterization of 
the s-wave deuteron wave function is 



1 N 



(41) 



.7 = 1 



with 



rn,j 



— a + (j — l)m . 

The usual boundary condition (finite wave function at the origin) leads to the constraint 

N 



0. 



(42) 



To put this realistic deuteron model on the light front, we work in momentum space and use the longitudinal 
momentum prescription above (Eq. ||) . The resulting form factor is completely similar to that of the Hulthen (Eq. 
pd| ) except there are now iV 2 terms instead of 4. At this point, the similarity leads us to suspect Q 4 behavior modified 
by a logarithm in asymptopia. Based on our above analysis, verification of the logarithm's presence requires that we 
check ip(x, p^) ~ (1 — x) 2 as x — > 1 and that the potential in r-space goes like 1/r near the origin. 

The momentum-space wave function has the end-point behavior 



N 



(l-x) 2 



; j (TO 2 + )' 



+ 0[(l~x) 3 } 7 



(43) 



where the term linear in 1 — x has vanished due to the constraint equation 42 , Appealing to the wave equation ^, 
we can determine the potential which generates this deuteron wave function 



Z)i=i C i 2 ( a ~ m o)m j + m 2 j 



V Mr (r) = 



2-2 



e — m„rj 



const, 



(44) 



where the constant ensures the potential vanishes as r — > oo. It is then straight forward to find the behavior near 
the origin 



lim V 



MT 



(r) 



2(a — TO ) + TO 



£(N) 



k(N) 



(45) 



where £(N) = ^2jCj and k(N) = ^j 2 Cj. Thus based on these analytic observations, equation |38| shows that even 
a realistic deuteron model will be troublesome in asymptopia. 



15 



7 Conquering Asymptopia 



As we have seen, equation £58] is a rather naive way to determine a form factor's asymptotic behavior. This equation 
and the analysis leading to it were extrapolated from our knowledge of non-relativistic wave functions. Indeed one 
may verify that equation (derived by an analysis parallel to the non-relativistic one PH ) gives exactly the same 
results as equation |38|. T he breakdown of factorization for the above models is clearly a relativistic problem (further 
verified by equation |37|). Equation ^8| may not be useful in determining the asymptotic behavior. Without it, 
however, we wouldn't be aware of the cause of our problems at high momentum transfer. Now knowing when to 
expect trouble in asymptopia, let us proceed to correctly deduce the asymptotic behavior. 



7.1 Wick-Cutkosky Model 

Before returning to the asymptotics of the Hulthen form factor, let us take a worthwhile look at the Wick-Cutkosky 
model. Consider two equally massive scalar particles which interact by exchanging a massless scalar particle. The 
potential for such a process has been found and consequently the ground-state wave function can be deduced using 
the momentum-space version of equation ^. The wave function is p7f 



ip(x,p A 



,5/2 



(P 2 



/^c 5 /V(l 



(4c(l - x)n 2 + p ±2 +{2x - l) 2 m 2 ) 2 (l + \2x - 1|) ' 



(46) 



where we have used the energy E{p) = \\J ^xii—™) * n ^ ne two-particle center of mass, and p 3 is given by equation 
|5|. Here k = \ma and a = g 2 /\6irm 2 with g as the coupling constant present in the interaction term. The invariant 
mass of the system is M = 2m — jma 2 . To write this wave function, we have converted the explicitly covariant 
form cited in p7{ into our own z-axis dependent form. The difference between these approaches does not concern 
us for the ground state of scalar particles. For a review of explicitly covariant light-front dynamics see p8| or p9[. 



The wave function in equation |46j is quite similar to our Coulomb wave function in section 3.1. The extra term in 
the denominator originates from retardation effects contained in the relativistic potential (effects which our naive 
models clearly lack) . The behavior of the wave function at the end point is not modified (to leading order in 1 — x) 
by retardation. Furthermore, the retarded potential is |2?| 

V(x,p x ;y,k ± ) WC = -47ra//C 2 , (47) 

with 



JC 2 = (p - k) 2 - (2.x - l)(2y - l)(E(p) - E(k)) 2 + 2\x - y\ (i?(p) 2 + E(ky 



M 2 



where (p — k) is given by equation 36. Using this potential for asymptotic Q, we note 



ni , (1 - I ) q ^V2.0)WC__^£_ (48) 

which is singular at x = 1. Given this and the wave function's end-point behavior, we once again appeal to equation 
pq and a logarithmically divergent x-integral confronts us in deducing the asymptotic behavior of the form factor. 
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Our experience above leads us to believe the true asymptotic behavior is Q~ A modified by a logarithm. This 
asymptotic behavior for the Wick-Cutkosky model has been found previously by Karmanov and Smirnov Q by 
considering regions which dominate the x and p 4 - integrals of the form factor. The same asymptotic behavior of 
the Wick-Cutkosky model was also found |)j by using the Bethe-Salpeter approach. Karmanov and Smirnov state 
that the logarithmic logQ 2 /Q 4 behavior was also been found earlier in |50| — a paper which admits the possibility of 
such logarithms only by announcing it is not considering such cases. The presence of logarithmic modifications to 
relativistic form factors is discussed in |l7| , where the authors interpret the Drell-Yan-West relation as valid modulo 
logarithms. Nonetheless, the correct asymptotic behavior of the Wick-Cutkosky form factor was deduced in [[)) as 
we shall now demonstrate using techniques considered above. 

Considering the analytic form of the wave function we can see a region which dominates the form factor for 
asymptotic Q: x near i. In this case, we can surely say (1 — x)Q « ^ ^> m and consequently the analysis leading to 
equation ^ is certainly valid. Moreover, the Wick-Cutkosky wave function is identical to our Coulomb wave function 
for x ~ ^. Thus appealing to equation we find (to leading order about x = i) 



lim F(Q 2 ) 

Q— >oo 



dx 



16m a 5 

ttQ 4 J Q f + (2x-l) 2 
16m 4 a 4 



[l + 0[a]), 



where we have assumed a <C 1 so that x(l — x)a 2 ~ a 2 /A. In order to compare with 
adopt their normalization (we have multiplied our expression ^|by ^(27r) -3 .) 



(49) 



we have been careful to 



The other dominant contribution in asymptopia comes from near the end-point region ^ < (1 



as) <C 1- A§ 



as we have seen above by producing logarithms from regularized x-integrals. We can thus deduce the remaining 
contribution in asymptopia via regularization. This result will be different from the Coulomb result, equation^, due 
to the retardation factor. First we note that in the near end-point region x(l — x)a 2 ~ 0. Now take 
and hence the contribution which interests us reads 



<x a < 1-AS 



lim F X (Q 2 ) = 



64m 4 a 5 
ttQ 4 

16m 4 a 5 
ttQ 4 



l-A 



x dx 



4(1 -x){2x- l) 2 



£(x , a, m) — log A + loj 



Q 



(l + O 



(50) 



Combining these two results (using only the logarithmic part of the latter), we arrive at the asymptotic behavior (to 
leading order in a) 



hm F(Q 2 ) w — - log^ + - 

Q->oo ttQ^ V m a/ 



(51) 



which agrees with the result found by Karmanov and Smirnov || . Furthermore one can use the Wick-Cutkosky wave 
function above to numerically calculate the form factor and test equation |5l| as predicting its asymptotic behavior. 
This form factor is less complicated than the Hulthcn model's. Consequently, the numerical integration is precise to 
larger Q. As before, we have graphically determined the coefficients a, (3 in equation ^9] and observed the asymptotic 
behavior as tending toward (a + /31og z)/z i , with z = Q. The coefficient f3 agrees well with equation ^l], differing 
by < 0.1% for a = 0.08,0.007. The error of the coefficient a depends on how rapidly the series in a converges. For 
a = 0.08, our graphically determined a differs from equation |l]by ~ 12%. While for a — 0.007, the error is ~ 2%. 
Indeed we have deduced the form factor's asymptotic behavior. 
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7.2 Back to the Hulthen 



Our analysis above has been quite general and we shall now apply it to the Hulthen form factor. To finish our quest 
through asymptopia, it remains to determine the coefficient a in equation 29 for the Hulthen model. As we learned 
above[], considering the contribution for x « | gives us a for asymptotic Q. So we return to equation ^3] and expand 
to leading order about x — h 



lim F(Q 2 



Anm 2 N{b 2 - a 2 ) 

Q 1 



log 



b 2 + (2x - l) 2 m 2 
a 2 + (2x - l) 2 m 2 



dx. 



(52) 



having used -C 1. At first glance, it appears we have dropped a factor of two from the asymptotic expression 

p3| . However careful consideration of region i (in equation ^l|) , shows its contribution vanishes (to leading order in 
1/Q). Evaluating the above integral and combining with the logarithmic part of our previous result (equation 28), 
we find 



lim F(Q 



2 _ 32irN(b 2 -a 2 ) 2 



Q 1 / rmr 

log — + o r— 

to 8 \ o + a 



(53) 



From which we deduce a — 0.046493 (and (3 = 0.045020 as discussed previously). Comparing with the graphically 
determined result of section 5.2, we see that there is a 17.8% difference. Again this difference is due to the series 
expansion in small parameters: a/m = 0.048943 and b/m = 0.29007. For smaller values of the parameters, we expect 
better results. However, with smaller parameters one needs higher Q to graphically determine a, (3 and the numerical 
integration becomes imprecise. Nonetheless, within our constraints we have verified equation 53 as the asymptotic 
behavior of the Hulthen form factor. 



8 Concluding Remarks 

We have undertaken a relatively simple task to compare relativistic and non-relativistic form factors for the Hulthen 
model of the deuteron. For small momentum transfer, the two versions differ by about a percent and the root mean 
square radii differ by even less. The behavior for large Q, however, led us on an unexpected journey. 

Our expedition through asymptopia helped us learn the Hulthen form factor's true behavior ~ logQ 2 /(5 4 for 
large Q. The path was circuitous because the conventional means (asymptotic expressions, Taylor series expansions) 
lead directly to logarithmic divergences. These difficulties are manifestations of the non-commutativity of the limits 
Q — ► oo and x — > 1, and hence indicative of the breakdown of factorization. 

Indeed, we find that this behavior is not particular to the Hulthen model. Equation ^38] tracks down the root 
of these difficulties. Generating our light front wave functions from non-relativistic potentials singular at the origin 
will lead to problematic relativistic form factors in the asymptotic limit. Such problems do not plague calculations 
in fundamental theories because higher-order corrections necessarily cancel the divergences. For realistic models, 
however, the breakdown of factorization persists and is an obstruction to straightforward asymptotic calculations. 
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